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ABSTRACT 

We present deep HST /ACS images and Keck spectroscopy of MC2 1635+119, a QSO hosted 
by a galaxy previously classified as an undisturbed elliptical. Our new images reveal dramatic 
shell structure indicative of a merger event in the relatively recent past. The brightest shells in 
the central regions of the host are distributed alternately in radius, with at least two distinct 
shells on one side of the nucleus and three on the other, out to a distance of '--^13 kpc. The light 
within the five shells comprises ^6% of the total galaxy light. Lower surface brightness ripples or 
tails and other debris extend out to a distance of ~65 kpc. A simple N-body model for a merger 
reproduces the inner shell structure and gives an estimate for the age of the merger between 
~30 Myr and ^^1.7 Gyr, depending on a range of reasonable assumptions. While the inner 
shell structure is suggestive of a minor merger, the total light contribution from the shells and 
extended structures are more indicative of a major merger. The spectrum of the host galaxy is 
dominated by a population of intermediate age (~1.4 Gyr), indicating a strong starburst episode 
that may have occurred at the time of the merger event. We speculate that the current QSO 
activity may have been triggered in the recent past by either a minor merger, or by debris from 
an older (~Gyr) major merger that is currently "raining" back into the central regions of the 
merger remnant. 



Subject headings: galaxies: active - galaxies: interactions 
quasars: individual (MC2 1635+119) 



galaxies: evolution — quasars: general 



^ Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated 
with program # GO-10421. 

^Department of Physics and Astronomy, Univer- 
sity of CaUfornia, Riverside, CA 92521, USA; email: 
gabriela. canalizoOucr . edu 

^Institute of Geophysics and Planetary Physics, Uni- 
versity of California, Riverside, CA 92521, USA; email: 
nicola.bennerfSucr.edu, bruno.jungwiert@ucr.edu 

^Astronomical Institute, Academy of Sciences of the 
Czech RcpubUc, Bocm' II 1401, 141 31 Prague 4, Czech 
Republic 

^Institute for Astronomy, University of Hawaii, 2680 
Woodlawn Drive, Honolulu, HI 96822, USA; email: stock- 
ton@ifa.hawaii.edu 

^Carnegie Observatories, 813 Santa Barbara Street, 
Pasadena, CA 91101, USA; email: schweizer@ociw.edu 

^Spitzer Science Center, California Institute of 
Technology, Pasadena, CA 91125, USA; email: 



1. Introduction 

The nature of QSO host galaxies has been 
debated for over four decades. Although the 
terms of the debate have gradually evolved dur- 
ing this time, there has been some progress. We 
now know, for example, that the majority of lu- 
minous low-redshift QSOs, whether radio loud 
or radio quiet, reside in the centers of galaxies 
that have relaxed light distributions like ellipti- 
c als (e.g.. iDisnev et al.lll995t iBahcall et al.lll997 : 



Dunlop et all l2003l : iFlovd et al.l 120041 ). This re- 
sult ties in nicely with the strong correlation, 
determined from galaxies with inactive black 
holes, between supermassive b lack hole mass and 
spheroid velocity dispersion ( Ferrarese fc MerrlO 
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20001: iGebhardt eralll2000l ): QSOs occur in the 
sorts of galaxies known to have the most massive 
central black holes. 

At the present epoch, only a tiny fraction of 
galaxies with massive spheroids shows luminous 
QSO activity. The very steep evolution of QSO ac- 
tivity with redshift indicates that some additional 
ingredient besides the mere presence of a super- 
massive black hole is necessary to produce QSO 
activity, and that this ingredient was much more 
common in the early history of the Universe. It 
has often been speculated that the mechanism un- 
derlying this evolution is the sudden inflow of gas 
to the center brought about by strong interactions 
or mergers. There has long been a fair amount of 
circu mstantial evide nce to support this idea (see, 
e.g., Stockton 11999 . and references therein), yet 
such arguments are by no means conclusive. 

The debate about the nature of QSO host 
galaxies presently centers on the question of how 
significant tidal interactions are for QSOs gener- 
ally: Do most QSOs at the current epoch begin 
their lives as mergers, or do most QSOs simply 
occur in old ellipticals to which nothing very in- 
teresting has happened recently? 

We are conducting a coordinated study with 
Keck spectroscopy and Hubble Space Telescope 
{HST) imaging of classical QSO host galaxies to 
investigate whether such hosts are truly quies- 
cent ellipticals with ancient stellar populations, or 
whether they are the results of mergers in the more 
recent past and have assumed elliptical morpholo- 
gies only as a result of violent relaxation due to 
the mergers. 

Elliptical hosts formed through mergers would 
be expected to show fine structure indicative of 
past tidal interactions, such as shells and rip- 
ples. Studies of nearby merger remnants (e.g. , 



Schweizer k Seitzeil Il992t ISchweizer et^ Il990l ) 



indicate that such structure can in general be de- 
tected even a few Gyr after the last major merger 
event. 

To look for any potential fine structure, we re- 
cently obtained very deep HST Advanced Cam- 
era for Surveys (ACS) images in a pilot study of 
five classical QSO host galaxies. In this paper, we 
present results for the first object, MC2 1635-1-119. 
The remaining four objects will be discussed in a 
subsequent paper (Bennert et al., in preparation). 



The host galaxy of MC21635-fll9 (z = 0.146; 
1"~ 2520 pc for = 0.7, l^^attor = 0.3, and 
Hf) = 71 kms~^ Mpc ~^) was first described by 
Hutchings et al. I (|l988h as having "slightly ellipti- 
cal amorphous structure" with a luminosity pro- 
file that does not follow a simple exponential or 
r^/^ law. Se veral companions are seen in the op - 
tical images (jHutching s et al. '1 19881: lMalkanl ll984l) 
as well as in the IR (Dunlo p et al .1119931). with- 



out a ny clear signs of i i iterac tion (jHutchings et al 



19881 ). iMcLure et all (|l999l ) compare fits to the 
host galaxy using an exponential disk and a de 
Vaucouleurs spheroid model, and conclude that 
the host resembles more close ly an elliptic a l. Re- 
garding the stellar contents, iNolan et al.l (|2001l ) 
estimate an age of 12 Gyr for the dominant stel- 
lar population in the host galaxy from off-nuclear 
spectra. 

Thus, previous studies seem to indicate that the 
galaxy hosting MC2 1635-1-119 is an elliptical with 
an old stellar population. We now present new 
HST and Keck observations that are in stark con- 
trast with any such conclusions. 

2. Observations and Data Reduction 



Spectroscopic observations and their analysis 
are described in detail elsewhere (Canalizo & 
Stockton, in preparation). Briefly, we obtained 
a spectrum of the host galaxy of MC2 1635-f 119 
with a total exposure time of 1.5 hours using 
the Low-Resolu tion Imaging Spectrometer (LRIS; 
Oke et al.|[l995l ) on the Keck I telescope on 2002 



March 4. We used the 400 groove mm ^ grism 
blazed at 3400 A for the blue side (LRIS-B), and 
the 300 groove mm~^ grating blazed at 5000 A 
for the red side (LRIS-R), yielding dispersions 
of 1.09 A pixel"^ and 2.55 A pixel"^ respectively. 
The slit was 1" wide, projecting to '^l pixels on 
the UV- and blue-optimized CCD of LRIS-B and 
~5 pixels on the Tektronix 2048x2048 CCD of 
LRIS-R. The slit position angle (PA) was 57°, 
placed roughly along the semi-major axis of the 
host galaxy and going through the QSO nucleus. 
The object was observed near transit, so that the 
effects of differential atmospheric refraction were 
minimized. 

The host galaxy spectrum was reduced using 
standard procedures. A scaled version of the QSO 
spectrum was subtracted from that of the host 
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galaxy; the spectrum was scaled by measuring the 
amount of flux in broad lines in the spectrum of 
the host. The final spectrum corresponds roughly 
to a region 2—5" from the nucleus on either side 
of the QSO and has a signal-to-noise ratio '^20. 
The spectrum was then modeled by performing 
least-squares fits t o the data using preli r ninary 
Chariot fc Bruzuall (j20o3) and iMarastonI ^QM) 
population synthesis models as described in § [6l 
Both the models and the observed spectrum were 
rebinned to the same spectral resolution. 

Imaging observations were obtained using 
ACS/WFC onboard the HST with the broad V- 
band F606W filter (AA = 2342A; 1 pixel corre- 
sponds to 0.05"). We obtained five sets of dithered 
images, each with four subsets of 550-586 s expo- 
sures, yielding a total integration time of 11432 s. 

We re-calibrated the data manually, starting 
from the pipeline flat-fielded individual exposures 
to improve the bias subtraction, i.e., to correct the 
offset (of a few DNs) between the adjacent quad- 
rants that is still present in the final product of 
CALACS ('Pavlovskv eTal]l2005l). W e then used 
MultiDrizzle (Koekemoer et al Jlooi) to combine 
the individual images, using the default values, 
bits=8578, as well as a deltashift-file contain- 
ing the offsets between the images as determined 
from stars within the field-of-view (FOV). The fi- 
nal distortion-corrected image is shown in Fig. [1] 
where the host galaxy shows clear shell structure. 



m 



3. Image Processing 

To enhance and analyze any fine structure that 
might be present, we applied various methods such 
as un sharp masking, creating a so-called structure 
map ( Pogge &: Martinill2002r ). as well as subtract- 



ing a central point spread function (PSF) for the 
QSO and a host galaxy model making use of GAL- 



FIT (jPeng et al.H2002f) . All different approaches 



confirm the existence of distinctive shells in the 
host galaxy (Fig. [2]), and we discuss each of them 
in turn. 

To create an unsharp-masked image, we divided 
the final image / by the / convolved with a Gaus- 
sian function of c = 5 pixel (G) : 



unsharp 



/ 



The structure map was derived by dividing / 
by the PSF-smoothed image (/ (g) P) and then 
convolving this ratio with the transpose of the PSF 
(P*): 



structure 



/ 



This process enhances unresolved or slightly re- 
solved structures on the scale of the PSF by remov- 
i ng the smooth light di stribution on larger scales 
(IPogge fc Martinill2002f ). 
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Fig. 1. — ACS/WFC image of MC2 1635-1-119, shown at different scales. Fine structure consisting of shells, 
arcs, and other debris is clearly seen at small and large scales. The images have been Gaussian smoothed 
with a sigma of either 0.5 pixels (left and central panels) or 2 pixels (right panel). In this and the following 
figures, north is up and east is to the left. 
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Fig. 2. — Different methods used to detect fine structure in MC2 1635+119, as described in the text. Top 
left: an unsharp-masked image, /unsharp- Top middle: a structure map, /structure- Top right: a residual image 
using GALFIT, where the model used for the host galaxy consists of a de Vaucouleurs profile only. Bottom: 
a residual image using GALFIT, where the model used for the host galaxy consists of a de Vaucouleurs and 
a Sersic profile of index n~l. 
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A PSF image is needed for both the struc- 
ture map and for modeHng with GALFIT. There- 
fore, we created an artificial PSF from TinyTim 
(Version 6.3) at the same position as our ob- 
ject as well as a "real" PSF using a star on an 
ACS/WFC F606W image. This image was ob- 
tained by searching the HST archive for a suitable 
star at roughly the same chip position as the QSO 
and with a high signal-to-noise ratio (S/N). We 
found a star at < 30 pixels away from the position 
corresponding to the QSO with a S/N of 20000 
that was observed on 20 dithered images with a 
total exposure time of '^SlOOs (GO-9433, datasets 
j6mfl9* and j6mf21*). We processed these im- 
ages in the same manner as described above for 
our data. 

In order to minimize introducing additional 
noise into the PSF subtraction and convolution 
operations, we first eliminated a few faint objects 
surrounding the PSF. Then, depending on the 
data values compared to the standard deviation 
s of the surrounding sky, we modified the PSF 
image as follows: (1) for data values > 7s, we 
retained the unmodified PSF; (2) for data values 
between 3 s and 7s, we smoothed the image with a 
Gaussian kernel with a = 0.5 pixel; (3) for smaller 
data values, we smoothed with a Gaussian kernel 
with a ~ 2.0 pixel; finally, (4) for data values that 
were < s after this last operation, we replaced the 
value with 0. 

To probe the quality of the two different PSFs, 
we subtracted them from both saturated and un- 
saturated stars within our FOV using GALFIT: 
the real PSF star gave significantly better results 
than the TinyTim PSF. From this exercise, we 
also determined the central region with a radius 
of ^1'.'7 is strongly affected by the PSF subtrac- 
tion because the QSO nucleus was saturated; any 
structure seen within this region is likely an arti- 
fact. 

The best enhancement of the shell structure was 
obtained using GALFIT (Fig. [2]) , a 2-dimensional 
galaxy fitting program capable of fitting simul- 
taneously one or more objects in an imag e with 
differe n t model light di s tribut ions (such as ISersic 
(1X968 ). de Vaucouleursl ( 19481 ). exponential, etc.; 
Peng et al.ll2002f) . Briefly, our adopted procedure 
was as follows: First, we created a mask to exclude 
the saturated pixels in the center, the diffraction 
spikes, any surrounding bright objects, and the 




shells themselves, in order to fit only the smooth 
underlying host galaxy light distribution. Then, 
a ("real") PSF as well as several Scrsic functions 
were fitted. In GALFIT, the Sersic power law is 
defined as 



S(r) = Eg exp 



where is the p ixel surface brightness at the ef- 
fective radius (jPeng et al. I [2002), and n is the 
Sersic index (n = 4 de Vaucouleurs, n = 1 expo- 
nential profile). In addition, we fitted the bright 
neighbor to the south of the QSO with a Sersic 
function. In all steps, the background sky was 
fitted simultaneously. This least-squares fit was 
then subtracted from the original image to gain 
the residual image, enhancing all structure that 
lies on top of the smooth host galaxy light distri- 
bution. 

When we used a single component for the host 
galaxy, the best fit was achieved with a Sersic func- 
tion of index n — 8.8. This fit was marginally bet- 
ter (only a few percent in y^) than the fit achieved 
using a de Vaucouleurs profile. On the other hand, 
the fit resulting from an exponential profile was 
much worse (roughly 50% in x^). This finding 
i s in a greement with the results bv lMcLure et al.l 
(Il999[ ). who determined that the host galaxy of 
MC2 1635-^119 is better fit by a de Vaucouleurs 
than an exponential profile. 

The fit improved substantially, however, when 
two components were included instead of one. 
Using two Sersic functions, the best result was 
achieved when one had an index n=4, which corre- 
sponds to a de Vaucouleurs profile, and the other 
an index n=0.91, which corresponds nearly to an 
exponential disk; this fit is shown in Fig. [2] and 
listed in Table [TJ If, instead, the index of one 
of the Sersic components was fixed to n=l (expo- 
nential), the best fit was achieved when the other 
component was close to a de Vaucouleurs profile 
(with index n=4.6; Table [J). Therefore, we con- 
clude that the host galaxy is well modeled by a 
de Vaucouleurs spheroid plus an exponential disk 
that makes up roughly one fourth of the light in 
the surface brightness profile, as detailed in Ta- 
ble [TJ In that table, we also list results for the 
fit using a de Vaucouleurs profile o nly in order to 



compa re our results with those of iDunlop et al 



(|2003l ). and we find that our results are very sim- 
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ilar to theirs. However, as Fig. [2] shows, the re- 
sulting model-subtracted image using only a de 
Vaucouleurs profile has residuals that are signifi- 
cantly larger than those obtained when we use a 
two-component model. 

4. Shell Structure and Luminosity 

In Fig. [3] we show a residual image of MC2 
1639-1-119 indicating the position of the different 
tidal features that we identify. The central circle 
with a radius of 1'.'7, corresponds to the area most 
affected by the saturated PSF; any features within 
this area may be artifacts of the PSF subtraction. 
Unfortunately, this prevents us from reliably de- 
tecting any shells or other structure that may be 
present in that region. 

The arcs labeled a through e in Fig. [3] are all 
segments of circles centered on the galaxy, empha- 
sizing the regularity of the interleaved shells. The 
projected radii of these shells are roughly 6.6, 7.6, 
8.3, 10.0, and 12.5 kpc, respectively. This set of 
bright shells is closely aligned with the semi-major 
axis of the host galaxy, at PA ^54°. The shell sys- 
tem shows roughly a biconical structure, although 
the edges of this putative bicone do not intersect 

HC2-1 635+1 19 

m 




5 arcsec 



at the center of the host galaxy. Shell e shows a 
discontinuity west of the QSO that may be due to 
obscuration by dust, or the shell may be made of 
two or more components. 

A set of lower surface-brightness shells or rip- 
ples (/, g and h) with seemingly different (greater) 
ellipticities are seen roughly perpendicular to the 
first set, both north-west and south-east of the nu- 
cleus. 

Further out to the north-east, there is an arc- 
like feature (i) extending out to a projected dis- 
tance of ~32 kpc. Other faint tails or wisps are 
seen in that same region (j). Finally, a much 
larger, faint and diffuse feature resembling either a 
shell or some tidal tail (fc) is visible ^65 kpc west 
of the nucleus. While this feature is very faint, 
we are confident it is real, particularly as this fea- 
ture is also visible in a WFPC2 archival image 
(GO-6776) when the image is median filtered and 
Gaussian smoothed. 

To estimate the luminosity within the shells 
compared to the total luminosity of the galaxy, we 
created a mask that includes all the light within 
an annulus of inner radius of 1'.'7 and outer radius 
of 6'.'8, but that at the same time excludes the 




* :■ 15 arcsec ^ 



Fig. 3. — Model-subtracted images of MC2 16354-119, where the most prominent fine-structure features are 
labeled. 
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Table 1 

Results of modeling the QSO host galaxy using GALFIT 



fit type 


function 


A(a,(5) (arcsec) 


nipeoew (mag) 


T'e (kpc) 


Sersic index 


hja 


PA (deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


de Vauc+Sersic 


SI 


(-0.03,0.04) 


17.46 


2.74 


4 (fixed) 


0.74 


57.3 




S2 


(0.72,0.26) 


18.80 


15.89 


0.91 (free) 


0.79 


28.8 


Sersic+Exp 


SI 


(-0.03,0.04) 


17.39 


2.75 


4.6 (free) 


0.75 


57.3 




S2 


(0.85,0.29) 


18.83 


16.5 


1 (fixed) 


0.79 


26.9 


de Vauc+Exp 


SI 


(-0.03,0.04) 


17.45 


2.68 


4 (fixed) 


0.74 


57.4 




S2 


(0.72,0.26) 


18.71 


16.03 


1 (fixed) 


0.8 


29.3 


de Vauc only 


SI 


(-0.03,0.04) 


17.26 


5.74 


4 (fixed) 


0.75 


52.3 


DunloD et al. (2003) 


SI 






5.73 


4 (fixed) 


0.69 


56 



Note. — Column (1) lists the GALFIT model, (2) the individual components used (S = Sersic), (3) the offsets with 
respect to the PSF, (4) the integrated apparent magnitude in the F606W filter, (5) the effective r adius, (6) the Sersic 
index, (7) the axis ratio, and (8) the position a ngle (east of north). Results from lPunlop et al.l | |2003D are listed for 
comparison. Note that the PA given here for the iDunlop et al.l (|2003| ) results was derived by adding the orientation of 
the spacecraft to the PA given in their Table 3, which was apparently not corrected for this orientation. 



diffraction spikes as well as several additional light 
sources from apparent companions. Note that the 
outer radius was chosen to be 3 x red with r^s 
determined from a single de Vaucouleurs fit (see 
Table [1]). This mask (good=l, bad=0) was mul- 
tiplied by the image and the total counts in the 
product were summed. This was done for both 
the GALFIT residual image (/shells) obtained sub- 
tracting the GALFIT model of a de Vaucouleurs 
+ exponential profile, and the GALFIT model 
itself (/galaxy)- Finally, we computed the ratio 
/shells //galaxy This yiclds the fractional luminos- 
ity of the shells between 1'.'7 and 6'.'8 radius as ~ 
6% of the host galaxy light (within the same an- 
nulus) . This estimated percentage may be smaller 
or larger than the true percentage depending on 
whether there are any shells within the radius af- 
fected by the PSF subtraction or not. 

Note that the percentage given refers to the to- 
tal flux within the shells out to 13 kpc. How- 
ever, the local contrast between the shells and the 
galaxy (as estimated by dividing the residual im- 
age by the GALFIT model) varies between 5 and 
20% and reaches 50% in Shell e. 

5. Time Constraints from Tidal Structure 

As described above, the host galaxy of MC2 
1635-1-119 reveals spectacular structure of regu- 
lar and aligned shells on projected radii of 5-13 



kpc. Similar she lls are observed in some local giant 



ellipt i cals (e.g., Malin &: Carter 1983t Schweizer 



IQSOl : Sikkema et al. 2007t ) andare interpreted as 



remnants of a merger event. It has been shown 
that the mergers tha t produce she ll-like structure 



can be either minor (lQuinnlll984l. hereaf ter Q84) 



or major (e.g.. lHernauist fc Spergellll992l ). In this 
section, we discuss both scenarios in the context 
of the morphology and physical size of the shells 
and structure we detect, with the aim of placing 
constraints on the age of the tidal interaction that 
formed them. 

5.1. Minor Merger 

We first consider the case of a minor merger 
since it allows for the simplest physical inter- 
pretation of the data. In this scenario, the 
system of regular concentric shells, confined 
within a finite range in azimuth, can result 
from the merger of a smaller galaxy (either spi- 
ral or elliptical) with a large elliptical along a 
nearly radial orbit ( Q84; Dupr az fc Combe"^ll986t 
Hernauist h Quinni [l988 . J98i hereafter HQ88 
and HQ89). 

The shell formation mechanism works as fol- 
lows: during the merger, stars from the smaller 
galaxy are captured by the massive galaxy and 
start to oscillate in its potential well. Since stars 
spend most of the time near the apocenters of their 
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orbits (where their radial velocities go to zero), 
a relative enhancement of the stellar density (a 
shell) forms there. The first shell is formed by 
captured stars that were initially in orbits with 
the smallest oscillation period, i.e., those with the 
smallest apocenter distance. 

As time goes on, the shortest-period stars move 
away from apocenter, while stars with slightly 
longer periods reach their apocenter at a slightly 
larger galactocentric distance. Due to a continu- 
ous range of oscillation periods, the first shell ap- 
pears to propagate radially outward while its stel- 
lar content progressively changes: it is thus a radi- 
ally propagating stellar density wave. A new trav- 
eling shell appears every time the shortest-period 
stars complete another oscillation period. After 
several oscillations, the massive elliptical galaxy 
reveals a system of shells where the outermost shell 
is the oldest, since this is the shell that formed 
first. This scenario gives a simple relation between 
the radius of this shell and the time of its forma- 
tion. 

We have constructed a simple N-body model 
that reproduces, at least qualitatively, the bright- 
est shells observed in MC2 1635+119. The N-body 
model uses the same technique as that used by Q84 
and HQ88. In this model, the secondary (smaller 
galaxy) moves on a radial orbit and is assumed to 
be disrupted instantaneously by the tidal forces 
of the primary (massive elliptical) after the first 
passage through the center of the primary. This 
corresponds to abruptly lowering the secondary's 
mass to zero, after which the test particles move in 
the potential of the primary alone. Thus, dynami- 
cal friction is assumed to be unimportant, and the 
model should only be considered as a zero-order 
description of the collision. 

We assumed a radial orbit with an initial sepa- 
ration between galaxies arbitrarily chosen to be 90 
kpc (5-18 times the scale-length of the primary). 
The initial infall velocity of the secondary was set 
equal to the escape velocity for the potential of 
the primary. 

We simulated the merger using (1) a de Vau- 
couleurs profile, and (2) a Plummer sphere (cor- 
responding to a Moffat's n=2 surface brightness 
profile). Since the goal of these simulations was to 
provide only a first order estimate of the merger 
timescale, we did not attempt to use more realis- 
tic composite density profiles of luminous and dark 



matter. We used effective radii ranging from 5 to 
20 kpc; this range s pans values for r^ff found by 



Dunlop et al. ( 2003) . ,Tavlor et al.l (|l996l ). and our 



own work (all corrected to the cosmology used in 
this paper). The mass of the giant elliptical was 
taken to be 3.2x10" Mq (jPunlop et al.ll2003l ). al- 
though we allowed for a range of masses up to 
3.2 X 10^^ M0 in order to account for a dark matter 
halo. The secondary-to-primary mass ratio and 
scale- length ratio were both fixed to 0.1; we note 
that while the precise choice of these two ratios 
is arbitrary, they affect mainly the contrast of the 
shells, and not the timescales, as long as the pri- 
mary dominates the potential. 

Figure 0] shows our results for simulations us- 
ing a Plummer surface brightness profile. The 
de Vaucouleurs model, which leads to lower con- 
trast and more spherical shells, will be discussed 
in more detail in a subsequent paper (Jungwiert 
et al., in preparation). Table [5] lists the timescales 
for two outermost shells (see below) to reach their 
observed radii in models with the range of parame- 
ters for the primary given above. We measure this 
timescale from the moment when the centers of 
mass of the two galaxies pass by each other (here- 
after "merger timescale"). We do not attempt 
to use the sizes or separations of inner shells to 
constrain the timescale since inner shells are more 
sensitive to the exact shape of the central density 
profile of the primary and are also more likely to 
be influenced by dynamical friction, which is not 
implemented in our model. 

Table 2 shows that, allowing for an uncertainty 
in the type of profile, for a rather large uncer- 
tainty in the effective radius, and for a consider- 
able amount of dark matter, the time for Shell e 
to reach its present distance of 12.5 kpc appears to 
be confined to a range of ^30—400 Myr after the 
centers of the two galaxies passed through each 
other. 

These ages are calculated assuming that Shell e 
is the outermost shell. However, we might consider 
the possibility that the tidal feature k may be a 
much older, fainter shell formed during the same 
encounter. This "shell" , at a projected distance of 
~65 kpc from the center of the host galaxy, would 
then give a merger timescale ranging from 100 Myr 
to 1.7 Gyr (see Table [2]), given the assumptions 
considered above. 

We emphasize that our simulations model the 
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Table 2 

Shell Formation Timescales from Numerical Simulations 



Rshell 


^^prirnary 


TdeV auc 


'^Plummer 


(kpc) 


(Mo) 


(Myr) 


(Myr) 


12.5 


3.2x10" 


100 - 245 


145 - 400 




3.2x10^2 


30 - 60 


45 - 135 


65 


3.2x10" 


360 - 1720 


1380 - 1620 




3.2x10^2 


100 - 400 


440 - 510 



Note. — The time range given for each model 
corresponds to a range of effective radii for the gi- 
ant eUiptical of 5 to 20 kpc. The time is measured 
from the moment when the centers of mass of the 
two galaxies pass through each other. 



1^ 



Fig. 4. — Shell structure in a restricted N-body 
simulation of a minor merger of two ellipticals 
(gE+dE). The masses of the galaxies are, respec- 
tively, 3.2 x 10" M© and 3.2 x 10^° M©, and their 
effective radii 5 and 0.5 kpc. Both galaxies are 
modeled as Plummer spheres. The smaller galaxy 
came from the right on a radial orbit. The box 
is 16x16 kpc. Only the particles belonging to the 
smaller galaxy are shown, to allow for comparison 
with images where a model of the host galaxy has 
been subtracted. 



simplest plausible case, and at this point we can- 
not exclude more complicated scenarios. In a sub- 
sequent paper (Jungwiert et al., in preparation) 
we will consider N-body simulations of this galaxy 
and of shell galaxies in general in more detail, fo- 
cusing on different gravitational potentials, var- 
ious mass ratios of colliding galaxies, dynamical 
friction, tidal stripping and the fate of gas. 

5.2. Major Merger 

While the numerical simulations described 
above can reproduce the morphology of the bright- 
est shells in MC2 1635-f 119, they do not rule out 
the possibility that the shells might have been 
created by a major merger. Further, the model- 
subtracted images (Figs. [5] and E]) show features 
(/i 5; h) that are off-axis from the direction of 
the encounter implied by the inner shells. Addi- 
tional tidal debris at different position angles is 
seen on much larger scales (features i, j, k). It 
is difficult to explain how all this structure might 
have formed as a result of a minor merger, pro- 
vided a single interaction is responsible for all the 
features. 

The fact that the inner shells appear to 
be closely aligned with the major axis of the 
host would also argue against a minor merger 

(see iHernguist fc Spergel 1992 . and references 

therei n). Using numerical simulations. lHernquist fc Spergel 
(jl992l ) show that mergers between two disk galax- 
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ies of similar mass can form shells, loops, and 
ripples. In particular, their simulations are 
compared to NGC 3923, one of the best exam- 
ples of a nearby ell iptical galaxy with shells 
IjMalin fc Carteil Il983l ). The system of shells of 
NGC 3923 (z = 0.005801) extends from distances 
close to the center (<2 kpc) out to ~100 kpc 
( Prieui 1988f ). The shells are distributed roughly 
in an hour-glass shape with an opening angle of ~ 
60°. While most of the shells appear aligned with 
the major axis of the galaxy, the outermost shell 
does not, a feature that is nicely reproduced by 
the simulations by Hernquist and Spergel. These 
characteristics are similar to those observed in 
MC2 1635+119, although it should be noted that 
the structure of the inner shells in MC2 1635+119 
is significantly more regular (non-intersecting and 
aligned) than that of the N GC 3923 shells or of 
the nu merical simulations bv lHernquist fc Spergell 
()l992ll . However, the comparison does point out 
that a major merger could also have formed the 
shells seen in MC2 1635+119. 

The amount of light observed in the shells may 
yield further clues to the nature of the merger. 
As mentioned in 21 the system of five bright 
shells comprises ~6% of the total luminosity of the 
galaxy. However, the shells contain only a fraction 
of the total number of stars that were originally 
part of the merging galaxy, i.e., those whose or- 
bital velocities are ne ar zero. Our numerical sim- 
ulations and those of iHernguist &: Spergell ( 19921 ) 
indicate that the stars in shells make up only one 
fourth or less of the total mass of the compan- 
ion. Therefore, assuming that the mass-to-light 
ratio is similar in both galaxies, the intruder may 
make up about 24% of the total mass. If we add 
to that the mass implied by the more extended 
"shell", the fraction may be closer to 30%. Thus, 
by this argument alone, the mass ratio of the origi- 
nal galaxies may have been close to 7:3 which may 
be considered a borderline major merger. 

Our simple N-body model produces shells up 
to a mass ratio of 3:1 for the parent galaxies. We 
did not investigate smaller mass ratios due to the 
increased complexity of such encounters. If we as- 
sume that the "shell" at 65 kpc (Shell fc) formed 
through a similar mechanism as that outlined in 
i)5.11 then the range of timescales of 100 Myr to 1.7 
Gyr would still hold for a major merger. If, how- 
ever, this feature was formed through the spatial 



wrapping of, e.g., a tidal tail, then estimating a 
timescale becomes more complex since timescales 
become more heavily dependent on initial condi- 
tions. As a reference, we note that simulat i ons o f 
the major merger in "The Mice" bv lBarne^ (2004) 
produce a merger remnant somewhat similar to 
MC2 1635+119 at a time close to 1 Gyr from the 
beginning of the merger event. 

6. Stellar Populations 

Figure [5] shows the Keck LRIS spectrum of the 
host galaxy of MC2 1635+119 in rest frame, rep- 
resenting its integrated light from 2" to 5" radius 
along the slit on either side of the nucleus (see [J2|) . 
Since the slit was placed roughly in the direction 
of the major axis of the host galaxy, the spectrum 
includes the brightest shells in the host (Fig. [3l a 
through e). 

The stellar component has a redshift z^hs = 
0.1474 (measured from absorption lines), equal 
to the redshift we measure from narrow emission 
lines, but slightly higher than that of the broad 
emission lines (z ~ 0.146). 

In order to model the sp ectrum, we used pop- 
ulation synthesis models by [Maraston' 1 20051) and 
the p reliminary models by Chariot fc Bruzual 
(l2007l) . We chose these two sets of models be- 




Fig. 5. — Keck LRIS spectrum of the host galaxy 
of MC2 1635+119 in rest frame. The black trace is 
th e observed spectruin . The red trace is the best 
fit Chariot fc Bruzual (l2007l ) model to the data. 
The model consists of 52% (by mass) of a 1.4 Gyr 
old population and 48% of a 12 Gyr population. In 
the bottom panel we show the residuals obtained 
by subtracting the model from the observed spec- 
trum. 
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cause they provide the best match to our spectral 
resolution and they both include contributions 
from thermally pulsating asymptotic giant branch 
stars, which are known to be particularly impor- 
tant in intermediate-ag e (~1 Gyr) stellar popu- 
lations (jMarastonl 120051 ) . Our original approach 
to analyzing the spectrum was to assume a domi- 
nant old stellar population (~12 Gyr) representing 
the population of the giant elliptical galaxy, with 
a smaller fraction of more recent star formation 
possibly triggered by the merger that formed the 
shells. Models that include a very small fraction 
(< 0.3%) of a young (< 50 Myr) starburst and a 
dominant ancient population can produce a rough 
fit to the continuum, but the fit to individual 
features such as Gall H&K and the CN band is 
rather poor. We tested spectral fits using differ- 
ent metallicities ranging from 0.02 to 2 solar, and 
found that solar metallicity models consistently 
yielded t he lowest y^. The choice of initial mass 



function ( Ghabr: 



St Y . 1 necnoice or in i tial mass 
ieJ booat iKroupal [200lt ISabeteil 



19551) made little or no difference 

However, the best fit to the observed spectrum, 
including both the continuum and stellar features, 
was achieved by adding a large contribution from 
an intermediate-age starburst population to the 
12 Gyr model. A better fit was achi eved with 
Chariot fc Bruzu"al ( 2007 ) than with iMaraston 
(|2005[ ) models, but both sets of models yielded 
simila r results. In the case of iGharlot fc Bruzuall 
(|2007l ) models, the best fit (shown in Fig. [5]) cor- 
responds to an intermediate-age population of 1.4 
Gyr contributing 52% of the tot al mass along the 
line of sight. The best fit using Marastoiil ( 20051 ) 
models is for an intermediate-age population of 
1.0 Gyr contributing 45% of the total mass along 
the line of sight. The real difference between the 
two models may be even smaller, considering that 
the Maraston models provide a coarser age grid 
(with steps in age at 1.0 and 1.5 Gyr) than the 
Chariot fc Bruzual models and the fact that the 
for the latter shows a shallow minimum from 
^1.2 to 1.9 Gyr (although the mass contribution 
from the starburst increases steeply with age). In 
both cases, increases rapidly beyond 2.0 Gyr. 

The determination of these intermediate-age 
components is robust with respect to the choice 
of the model for the older population: the same 
intermediate-age populations are obtained when 
the older population is varied from 6 to 14 Gyr. 



If we use models of metallicities lower than solar, 
a single population can be used to fit the data, 
although the overall fit is significantly worse. In 
this case, the oldest population that yields a rea- 
sonable fit is less than 3 Gyr old. Single popula- 
tions older than 4 Gyr yield poor fits regardless 
of their metallicity or initial mass function. Al- 
though it is possible that the spectrum may be 
somewhat reddened by dust, it is unlikely that the 
age of the starburst component would be signifi- 
cantly younger than one Gyr, given the absorption 
features that we observe. Finally, an inaccurate 
subtraction of the QSO contribution could affect 
the shape of the continuum. We tested the effects 
of this by fitting spectra that were slightly over- 
subtracted and under-subtracted. While the 
for these cases was somewhat larger, the age of 
the starburst for the best fit remained the same. 

Naturally, the number of possible combinations 
of populations to model the spectrum of the host 
is large. We have kept our analysis simple by test- 
ing only a limited number of possibilities corre- 
sponding to physically plausible scenarios. There- 
fore, while we cannot exclude more complex star- 
formation histories, we are fairly certain that: (1) 
The dominant component of the stellar population 
in the host of MG2 1635-f 119 is not ancient, and 
(2) A small percentage by mass of recent (less than 
a few hundred Myr) star formation superposed on 
an old (> 6 Gyr) population can be ruled out, re- 
gardless of the age of the dominant population. In- 
stead, the spectrum of the host of MC2 1635-1-119 
is dominated (at least in flux) by an intermediate 
age population of 1-2 Gyr. 

7. Discussion 

In agreement with previous observations, we 
have found that the surface brightness profile of 
the galaxy hosting MC2 1635-1-119 is closer to a de 
Vaucouleurs than an exponential profile. However, 
our new ACS image reveals that a fainter exponen- 
tial profile is also present, comprising up to one 
forth of the total luminosity. Moreover, our ob- 
servations have uncovered a spectacular system of 
shells and other faint structure in the host galaxy 
at small and large scales, showing that the host is 
far from being undisturbed. We have also found 
that the stellar populations in the host galaxy 
seem to have a substantial contribution (~50% 
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by mass) of an intermediate-age stellar population 
from a 1-2 Gyr old starburst. 

While the large contribution of an intermediate- 
age population to the spectrum of the host galaxy 
of MC2 1635-1-119 is intriguing, it is by no means 
unusual. Recent studies of AGN host galax- 
ies fe.g.. IJahnke et alJ [ 20041: ISanchez et all 120041 : 



Kauffmann et alJ l2003t ICanalizo et all l2006 fMn- 
dicate that galaxies hosting the most luminous 
AGN are often dominated by bulges whose col- 
ors are significantly bluer than those of inactive 
elliptical galaxies and are consistent with the pres- 
ence of intermediate-age starbursts. Based on po- 
sitions of the hosts in the D„(4000)/H5^ plane, 
iKauffmann et al. ( 20031) suggest that these AGN 
hosts have had significant bursts of star formation 
in the past 1-2 Gyr. 

Why do AGN host galaxies show these strong 
intermediate-age populations? And, what is the 
physical connection, if any, between the putative 
~Gyr old starburst and the nuclear activity? Un- 
derstanding the nature of this relation is impor- 
tant because it could have implications for the 
triggering mechanisms and duty cycles of AGN. 
Our study of MC2 1635-1-119 provides some clues 
that may be applicable to a larger population. 

We now know that the host galaxy of MC2 
1635-1-119 was unequivocally involved in a tidal 
encounter. Our rough estimates discussed in § [5] 
place the timescale for this encounter at less than 
~1.7 Gyr, which could be compatible with the age 
of the major starburst. However, the large uncer- 
tainty in our estimate does not rule out the possi- 
bility of a substantially more recent event. We are 
also unable to discriminate between a major and 
a minor merger as the culprit for the shell struc- 
ture that we observe. Our results give us enough 
information, however, to speculate on a couple of 
likely scenarios. 

First, consider the case where the inner shell 
structure was formed through the accretion of a 
low-mass companion (one tenth or less of the mass 
of the primary). The overall morphology that we 
observe would have to be caused by more than 
one event, and the fact that there was a dra- 
matic episode of star formation more than one 
Gyr ago would argue for a past (major?) merger 
connected to the large-scale tidal debris. In that 
case, it is possible that the giant elliptical pos- 
sessed a higher gas content as a result of the past 



merger event, and so the QSO activity was more 
readily triggered (or rejuvenated) in it by a mi- 
nor merger than it would have been in a gas- 
poor elliptical would. This may well be the case 
in CygnusA, where an ongoing minor merger ap- 
pears to be responsible for tr iggering the nuclear 
activity (jCanalizo et al. l2003l ). 

Consider now the alternative case where a ma- 
jor merger is responsible for both the starburst 
and all of the structure that we observe. This 
merger event would have occurred over one Gyr 
ago and would have likely (though not necessar- 
ily) triggered a first episode of accretion onto the 
black hole(s). Feedback from the QSO quenched 
any further star formation. Assuming theoretical 
estimates for the durati on of QSO activity are cor- 
rect (e.g., 10^ -10^ vr: lYu fc Tremaindl20oi ) . the 
QSO activity would have ceased as the merger con- 
tinued its course and the morphology of the newly 
merged galaxies began to relax into the shape of 
an elliptical. Eventually, the extended tidal debris 
would "rain" back into the central regions of the 
galaxy, triggering a new episode of QSO activity. 
A time delay in the onset of QSO activity would be 
in agreement with predictions by hydrody namical 
simulati ons of mergi ng galaxi es (see e.g.. iBarned 
i998; Springel et all 12005: Hopkins et all |2007[) . 
These models frequently predict a second peak in 
star formation that also occurs much later in the 
merger. Since our spectroscopic observations ex- 
clude a radius of ~5 kpc around the nucleus, we 
would not have detected any recent star formation 
that may be present in the central regions of the 
host galaxy. 

While these are interesting scenarios, they are, 
for the moment, no more than "guided" specula- 
tion. More complete N-body models as well as 
high angular-resolution spectroscopy to measure 
the kinematics of the stellar component are needed 
to get a better handle on the kind of encounter 
that formed the observed structure. However, we 
will also need to study larger samples to attempt 
to answer more complex questions, such as the 
precise timing of the triggering of the QSO activ- 
ity, which in turn should help answer questions 
regarding duty cycles and feedback. 
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